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The arrival time probability distribution is defined by analogy with the classical mechanics. The 
difficulty of requirement to have the values of non-commuting operators is circumvented using the 
concept of weak measurements. The proposed procedure is suitable to the free particles and to the 
particles subjected to an external potential, as well. It is shown that such an approach imposes an 
inherent limitation to the accuracy of the arrival time determination. 
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^ ■ I. INTRODUCTION 

Time plays a special role in quantum mechanics. Unlike other observables, time remains a classical variable. It 
is well-known, that the self-adjoint operator of time does not exist for bounded Hamiltonians. The problems with 
' the time rise from the fact that in quantum mechanics many quantities cannot have definite values simultaneously. 
^ , However, the detection of the particles in time-of-flight and coincidence experiments are common, and quantum 
>^ • mechanics should give the method for the calculation of the arrival time. The arrival time distribution may be useful 
' solving the tunneling time problem, as well. Therefore, the quantum description of the arrival time has attracted 
O ■ much attention 0. ^- ^ *4 6, 2, J. 1, IQ., 1 h 11 E lH d • 
' Aharonov and Bohm introduced the arrival time operator m 

: 
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' By imposing several conditions (normalization, positivity, minimum variance, and a symmetry with respect to the 
Oh, arrival point X) a quantum arrival time distribution for the free particle has been obtained by Kijowski 2]. The 
-j-J ' Kijowski's distribution may be associated with the positive operator valued measure generated by the eigenstates of 
^ , Tab- However, Kijowski's set of conditions cannot be a ppl ied in a general case 0. Nevertheless, the arrival time 
^ operators can be built even if the particle is not free [irI Ii3 | . 

In this paper we take another approach. Since the mean arrival time even in the classical mechanics can be infinite 
• • ■ or the particle may not arrive at all, it is convenient to deal not with the mean arrival time and corresponding operator 
T, but with the probability distribution of the arrival time. The probability distribution of the arrival time can be 
obtained from the suitable classical definition. The non-commutativity of the operators in quantum mechanics is 
circumvented by using the concept of weak measurements. Such an approach has several advantages. It gives, in 
principle, the procedure of measuring the obtained quantity. Since in the classical mechanics all quantities can have 
definite values simultaneously, weak measurements give the correct classical limit. The concept of weak measurements 
has been already applied to the time problem in quantum mechanics . 

We proceed as follows. In Sec. ^ we discuss various definitions of the arrival time in the classical mechanics. The 
weak measurement of the quantum arrival time distribution is presented in Sec. IIIII The properties of the obtained 
quantity are analyzed in Sec. II VI Section Ivl summarizes our findings. 
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II. ARRIVAL TIME IN CLASSICAL MECHANICS 



In classical mechanics the particle moves along the trajectory H{x,p) — const as t increases. This allows us to work 
out the time of arrival at the point x{t) = X , hy identifying the point (a;o,po) of the phase space where the particle 
is at i = 0, and then following the trajectory that passes by this point, up to the arrival at the point X. If multiple 
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crossings are possible, one may define a distribution of arrival times with contributions from all crossings, when no 
distinction is made between first, second and nth arrivals. In this article we will consider such a distribution. 

We can ask, whether there is a definition of the arrival time, which can be valid both in classical and quantum 
mechanics. In our opinion, the words "the particle arrives from the left at the point X at the time t" mean that: (i) 
at time t the particle has been in the region x < X and (ii) at time t + At {At 0) the particle has been found in 
the region x > X. Now we apply this definition, given by (i) and (ii), to the time of arrival in the classical case. 

Since quantum mechanics deals with probabilities, it is convenient to use probabilistic description of the classical 
mechanics, as well. Therefore, we will consider an ensemble of non interacting classical particles. The probability 
density in the phase space is p{x,p]t). 

Let us denote the region x < X dsTi and the region a; > X as The probability that the particle arrives from 
the region Fi to the region T2 at the time between t and t + At \s proportional to the probability that the particle is 
in the region Fi at time t and it is in the region F2 at time t + At. This probability is 

n+(i)At = -^ / dpdxpix,p;t), (2) 

where is the constant of normalization and the region of the phase space ft has the following properties: (i) the 
coordinates of the points in f2 are in the space region Fi and (ii) if the phase trajectory goes through the point of the 
region at time t then the particle at time t + At is in the space region F2 . Since At is infinitesimal, the change of 
the coordinate during the time interval At is equal to ^At. Therefore, the particle arrives from the region Fi to the 
region F2 only if the momentum of the particle in the point X is positive. The phase space region 51 consists from 
the points with positive momentum p and with coordinates between X ^ p/mAt and X. Then from Eq. (|2Jl we have 
the probability of arrival time 

Il+{t)At=— dp dxp{x,p;t). (3) 

Since At is infinitesimal and the momentum of every particle is finite, we can replace x in Eq. JSJ by X and obtain 
the equality 



U+it,X) = ^f dp^p{X,p;t). 
N+ Jo m 



(4) 



The obtained arrival time distribution Il+{t,X) is well-known and has appeared quite often in the literature (see, 
e.g., review |l5l | and references therein). 

The probability current in classical mechanics is 

/■+°° p 

J{x;t)^ / —p{x,P]t)dp. (5) 

From Eqs. (@J and (O it is clear that the time of arrival is related to the probability current. This relation, however, 
is not straightforward. We can introduce the "positive probability current" 

/■°° P 

J+{x;t)= —p(x,p;t)dp (6) 
Jo m 

and rewrite Eq. Q as 

n+{t,X)^^J+{X;t). (7) 

The proposed [11 El 111 various quantum versions of J+ even in the case of the free particle can be negative (the 
so-called backflow effect). Therefore, classical expression Q for the time of arrival becomes problematic in quantum 
mechanics. 

Similarly, for the arrival from the right we obtain the probability density 



n_(t,X) = ^J_(X;i), (8) 



where the negative probability current is 



J-{x;t)^[ —p{x,p;t)dp. (9) 
J-00 m 
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We see that our definition, given at the beginning of this section, leads to the proper result in the classical mechanics. 
The conditions (i) and (ii) does not involve the concept of the trajectories. We can try to use this definition also in 
quantum mechanics. 



III. WEAK MEASUREMENT 



The proposed definition of the arrival time probability distribution can be used in the quantum mechanics only if 
the determination of the region, in which the particle is, does not disturb the motion of the particle. This can be 
achieved using the weak measurements of Aharonov, Albert and Vaidman 

We have the detector in the initial state |$). The detector interacts with the particle only in the region Fi. In 
order the weak measurements can provide the meaningful information, the measurements must be performed on an 
ensemble of identical systems. Each system with its own detector is prepared in the same initial state. The readings 
of the detectors are collected and averaged. 

We take the operator describing the interaction between the particle and the detector of the form fvA 

Hi = XqPi, (10) 

where Pi is the projection operator, projecting into the region Fi and A characterizes the strength of the interaction 
with the detector. The interaction operator H10() only slightly differs from the one used by Aharonov, Albert and 
Vaidman 21] . The similar interaction operator has been considered by von Neumann and has been widely used 
in the strong measurement models (e.g., [s^, |^ 113, 113, 13 many others). 

The measurement time is r. We assume that the interaction strength A and the time t are small. A very small 
parameter A ensures the undisturbance of the particle's evolution. The operator q acts in the Hilbert space of the 
detector. We require a continuous spectrum of the operator q. For simplicity, we can consider this operator as the 
coordinate of the detector. The momentum, conjugate to q, is pq. 

Since interaction strength A and the duration of the measurement r are small, the probability (Pi) to find the 
particle in the region Fi does not change significantly during the measurement. The action of the Hamiltonian Hl()(l 
results in the small change of the mean detector's momentum (pq) — {pq)o — — At(Pi), where {pq)o = ($(0)|pg|<f>(0)) 
is the mean momentum of the detector at the beginning of the measurement and (pq) — (<I'(T)|pg|$(r)) is the mean 
momentum of the detector after the measurement. Therefore, in analogy to Ref. (2lf . we define the "weak value" of 
the probability to find the particle in the region Fi, 

w^(i) = (A) = -^^^l^—M. (11) 

Ar 

In order to obtain the arrival time probability using the definition from Scc.^ we measure the momenta Pq of each 
detector after the interaction with the particle. After the time At we perform the final, postselection measurement 
on the particles of our ensemble and measure if the particle is found in the region F2. Then we collect the outcomes 
Pq only for the particles found in the region F2. 

The projection operator projecting into the region F2 is P2. In the Heisenberg representation this operator is 

P2{t) - U{t)^P2U{t), (12) 

where tl is the evolution operator of the free particle. After the measurement the state of the particle and the 
detector is C/m(t)|$) Ivf), where is the initial state of the particle and J7m is the evolution operator of the particle 
interacting with the detector. The joint probability that the detector has the momentum pq and the particle after the 
time At is found in the region F2 is 

W{pq,2) - {^mUM{ry\Pq){Pq\P2{At)UM{r)m\^}, (13) 

where \pq) is the eigenfunction of the momentum operator j5g, P2 is the projection operator, projecting into the region 
F2. In quantum mechanics the probability that two quantities simultaneously have definite values does not always 
exist. If the joint probability does not exist then the concept of the conditional probability is meaningless. However, 
in our case operators P2 and \pq){pq\ act in different spaces and commute, therefore, the joint probability W{pq,2) 
exists. 

Let us define the conditional probability, i.e., the probability that the momentum of the detector is Pq, provided 
that the particle after the time At is found in the region F2. This probability is given according to the Bayes's theorem 



W{pq,2 ) 

W{2) 



mP,\2) - (14) 
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where W{2) is the probability that the particle after the time At is found in the region The average momentum 
of the detector on condition that the particle after the time At is found in the region r2 equals to 

{p,)2 - J p,W{p,\2)dp, = ^(*|($|;7M(r)tp,P2(Ai);7M(r)|$)|vI/), (15) 
where P2{At) — U{AtyP2U{At) is the operator P2 in the Heisenberg picture. According to Eq. (|ll(l 

|.^ ^ {P,)0-{P,)2 ^ {P,)0W{2) - {p,)2W{2) 

^ ' ^ Ar AtVK(2) ^ ' 

is the weak value of probability to find the particle in the region Fi on condition that the particle after the time At is 
in the region r2 . The probability that the particle is in the region Fi and after the time At is in the region F2 equals 
to 

W{1,2) = W{2)W{1\2). (17) 

When the measurement time t is sufficiently small, the influence of the Hamiltonian of the particle can be neglected 
and the evolution operator of the particle and the detector can be expressed as 

z 

Um{t) = exp ( --XrqPi 

We expand the operator Um{t) into the series of the parameter A, assuming that A is small. In the first-order 
approximation, using Eqs. (|16|l and H17|) . we obtain 

W{1,2) « i(P2(A0A + AP2(At)) + {{pg){q)-Re{qpg)) ([A, A(Ai)]). (18) 

The probability 1^(1,2) is constructed using the conditions (i) and (ii) from Sec. ^1 the weak measurement is 
performed to determine if the particle is in the region Fi and after the time At the strong measurement determines 
if the particle is in the region F2. Therefore, according to Sec.m the quantity W{1,2) after normalization can be 
considered as the weak value of the arrival time probability distribution. 

Eq. ({TH|l consists of two terms and we can introduce two quantities 

n^'^ = ^^^1^2 (AO + A (AO A) (19) 



and 



Then 



n('^-^([A,A(At)]). (20) 



9 At 

w^(i,2) -n(i)Ai- ^((A)(<7) -RG(#,))n(2). (21) 

If the commutator [Pi, A (At)] in Eqs. (|19|I - H21|) is not zero then, even in the limit of the very weak measurement, 
the measured value depends on the particular detector. This fact means that in such a case we cannot obtain the 
definite value for the arrival time probability. Moreover, the coefficient ({pq){q) —Iie{qpq)) may be zero for the specific 
initial state of the detector, e.g., for the Gaussian distribution of the coordinate q and momentum pq. 

The quantities 1^(1,2), n(i) and H^^) are real. However, it is convenient to consider the complex quantity 

He = n(i) + zn(2) = i-(AA(Ai)). (22) 

We name it "complex arrival probability" . We can introduce the corresponding operator 

n+ = ^AA(AO. (23) 



5 



By analogy, the operator 

n_ = l-^P,P,{At). (24) 

corresponds to the arrival from the right. 

The introduced operator 11+ has some properties of the classical positive probability current. From the conditions 
A + A = 1 and Pi{At) + ^(Af) = 1 we have 

n+-n_ = i^(A(At)-A). 

In the limit At ^ we obtain the probability current, J = limAt^o(n+ — 11-), as it is in the classical mechanics. 
However, the quantity (n_|_) is complex and the real part can be negative, in contrary to the classical quantity J+. 
The reason for this is the non-commutativity of the operators Pi and A (At). When the imaginary part is small, the 
quantity (n_|_) after normalization can be considered as the approximate probability distribution of the arrival time . 



IV. ARRIVAL TIME PROBABILITY 

The operator 11+ is obtained without specification of the Hamiltonian of the particle and is suitable to the free 
particles and to the particles subjected to an external potential as well. In this section we consider the arrival time 
of the free particle. 

The calculation of the "weak arrival time distribution" W{1,2) involves the average (11+) . Therefore, it is useful 
to have the matrix elements of the operator 11+ . It should be noted, that matrix elements of the operator 11+ as well 
as the operator itself are only auxiliary and does not have independent meaning. 

In the basis of momentum eigenstates normalized according the condition {pi\p2) — 2'KhS{pi — P2), the matrix 
elements of the operator 11+ are 

(Piin+b) = ^(pi|A?7(At)tA?7(At)|p2) 

= — / dxi da;2e-^fi^i(a;i|;7(At)^|a;2)e^f^^=-^^^*. (25) 

At J_oo Jx 

After performing the integration we obtain 



(Pi|n+b2) = TTTTT— rexp ( ^(p2 -Pi)^ 



where y/i = exp(i7r/4). When 



1 At , o At At 



the matrix elements of the operator 11+ are 

(pi|n+b2)«^i^exp(^J(p2-Pi)x). (27) 

This equation coincides with the expression for the matrix elements of the probability current operator. 
From Eq. (jSEJ we obtain the diagonal matrix elements of the operator 11+ 



{p\n+\p) = ^ erfc ( ) + , ^ e-^^^*. (28) 
^ ' ' 2m \^ 2hm J ^t2nhmAt ^ ' 

The real part of the quantity (p|n+|p) is shown in Fig. ^ and the imaginary part in Fig.|21 
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FIG. 1: The real part of the quantity {p|li+|p), according to Eq. 1281 . The corresponding classical positive probability current 
is shown with the dashed line. The used parameters are h — 1, m = 1, and At — 1. In this system of units, the momentum p 
is dimensionless. 

Using the asymptotic expressions of the function crfc wc obtain from Eq. H28|l that 

lim {p\n+\p)^P- 
p— (+00 m 

and (p|n+ \p) — !■ 0, when p —^ — oo, i.e., the imaginary part tends to zero and the real part approaches the corresponding 
classical value as the modulus of the momentum \p\ increases. Such a behaviour is evident from Figs. ^ and |21 too. 

The asymptotic expressions of the function erfc are valid when the argument of the erfc is large, i.e., \p\\J^^^ > 1 

or 



A.>A. (29) 

Here is the kinetic energy of the particle. The dependence of the quantity Re(p|n_|_||5) from At is shown in Fig. 
For small At the quantity (p|n+|p) is proportional to l/y/At. Therefore, unlike the classical mechanics, in quantum 
mechanics At cannot be zero. Eq. (|29|l imposes the lower bound to the resolution time At. It follows that our model 
does not permit determination of the arrival time with the resolution greater than h/Ek- The relation, similar to 
Eq. H29|l . based on measurement models has been obtained by Aharonov et al. |29|. The time-energy uncertainty 
realtions, associated with time of arrival distribution, are also discussed in Refs. I^, |33 



V. CONCLUSION 



The definition of density of one sided arrivals is proposed. This definition is extended into quantum mechanics, 
using the concept of weak measurements by Aharonov et al. The proposed procedure is suitable to the free particles 
and to the particles subjected to an external potential, as well. It gives not only the mathematical expression for 
the arrival time probability distribution but also the way to measure the obtained quantity. However, this procedure 
gives no unique expression for the arrival time probability distribution. 

In analogy with complex tunneling time the complex arrival time "probability distribution" is introduced (Eq. I|22|) '). 
It is shown that the proposed approach imposes an inherent limitation, Eq. (|29|l . to the resolution time of the arrival 
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FIG. 2; The imaginary part of the quantity (pjri+|p). Used parameters are the same as in Fig. 
time determination. 
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FIG. 3: The dependence of the quantity Re(p|li+|p) according to Eq. (I28I I from the resolution time At. 
classical positive probability current is shown with the dashed line. The used parameters are h = 1, m — 1, 
units, the time At is dimensionless. 
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